Astronomy & Astrophysics manuscript no. clref.tex 


February 5, 2008 


(DOI: will be inserted by hand later) 





ISM gas removal from starburst galaxies and the premature 

death of star clusters 

Claudio Melioli^ and Elisabete M. de Gouveia Dal Pino^ 



in 
o 
o 

> 
O 

:z; 

(N 



> 
m 



in 
o 

Oh 

6 



Universidade de Sao Paulo, TAG, Rua do Matao 1226, Cidade Universitaria, Sao Paulo 05508-900, Brazil 
e-mail: cmelioli@astro.iag.usp.br , dalpino@astro.iag.usp.br 

Received ; accepted 

Abstract. Recent observational studies of the age distribution of star clusters in nearby merging galaxies and 
starburst (SB) galaxies indicate a premature death of the young clusters. The fate of an evolving star cluster 
crucially depends of its gas content. This behaves like a glue that helps to keep the star system gravitationally 
bound. In SB systems where the rate of supernovae (SNe) explosions is elevated one should expect an efficient 
heating of the gas and its complete removal which could then favor the rapid dissociation of the evolving star 
clusters. Based on a contemporaneous study of the dynamical evolution of the interstellar gas in SB environments 
(Melioli & de Gouveia Dal Pino 2004) where it has been considered also the presence of dense clouds that may 
inhibit the heating efficiency of the interstellar gas by the SNe, we have here computed the timescales for gas 
removal from young clusters embedded in these systems and found that they are consistent with the very short 
timescales for cluster dissolution which are inferred from the observational studies above. Our results indicate 
that typical SB proto-clusters should start to disperse after less than 5 Myr. For a given total gas mass content, 
this result is nearly insensitive to the initial star formation efficiency. 
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1. Introduction 

Star clusters in nearby SB galaxies seem to contribute 
with only ~20% of the light seen in UV images, so that 
the diffuse light is dominated by the contribution of field 
stars (Meurer et al. 1995). There is increasing observa- 
tional evidence that suggests that this diffuse UV field 
light is actually due to the dispersion of aged star clusters 
that have left their remaining stars in the field (Tremonti 
et al. 2001, Lada & Lada 2003, Fall 2004). Chandar et 
al. (2005) have observed this trend in eleven nearby SB 
galaxies and have also identified a lack of massive stars 
in the field which supports the interpretation above that 
the field stars are composed of older, dissolving clusters. 
They also estimate that the star clusters need to dissolve 
on timescales 7-10 Myr in order to create the field stars. 

A contemporaneous statistical study of the age dis- 
tribution of star clusters in nearby merging galaxies (the 
Antannae galaxies) by Fall (2004; see also Whitmore et 
al. 2005) has shown that the number of clusters decreases 
by a factor of ^ 10 by the time the cluster population has 
reached the age of 10 Myr. This rapid decline is very 
insensitive to the cluster mass, at least for masses above 3 
xlO"* solar masses, and should be an indication that most 
of the clusters have become gravitationally unbounded. 
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The fraction of gas mass that is converted into stars 
in SBs is about 10 times larger than in normal galaxies 
and, as such, the production of a high number of mas- 
sive stars results high rates of SN explosions. Thus SN 
explosions, as well as the ionizing radiation from massive 
stars, and stellar winds may efficiently heat the gas and 
eventually remove a significant fraction of the ISM from 
a proto-cluster, leaving the stars gravitationally unbound. 
The crossing time for stars orbiting within a bound clus- 
ter or proto-cluster is Tcr — 2R/a, where a is the velocity 
dispersion of cluster, a ~ [OAG{Mstar + Mgas) / R]^^"^ , and 
R is its the radius. Now, if a proto-cluster suddenly loses 
most of its mass by the removal of the gas and the in- 
stantaneous velocity dispersion of the stars at the time 
of gas removal, a, becomes larger than the escape veloc- 
ity, Vesc — [{'2GAdstar) / R]^^"^ , then the cluster will be no 
longer gravitationally bound and will expand almost freely 
with its radius increasing with age as (see, e.g.. Fall 2004; 
Lada & Lada 2003) Rt — R{T/Tcr) and its surface density 
decreasing as ~ Hir/Tcr)- 

A recent independent theoretical study of the evolution 
of the ISM of SB environments by Melioli & de Gouveia 
Dal Pino (2004; see also Melioli, de Gouveia Dal Pino & 
Raga 2005), has revealed that the effectiveness of the SNe 
and their shock waves to heat the IS gas and make it to 
be ejected from the SB system is sensitive to the total gas 
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mass content. In particular, we have found that if the total 
gas mass both in the form of clouds and diffuse gas is large 
enough, then the radiative cooling time scale of the gas 
can remain shorter than the timescale for the development 
of a hot superbubble (due to the SNRs interactions) for 
about half lifetime of the SB. During this time, most of 
the gas is retained in the SB, but after this period the 
SNe heating efficiency rapidly increases to its maximum 
value thus leading to gas expansion and removal from the 
SB in a short period. The clouds play a fundamental role 
in this process as they work like valves that are able to 
both retain part of the gas to themselves and lose part of 
it to the diffuse ISM through cloud photoevaporation and 
then, maintain the ISM density and the radiative cooling 
rate high enough for a long time. As long as the gas is 
retained it may favor new generations of star formation. 

In this work, we apply the analysis above, which was 
performed for an entire SB system (Melioli & de Gouveia 
Dal Pino2004), to the scale of individual star proto- 
clusters embedded in it, and show that the resulting typ- 
ical timescales both for gas retention and removal from 
these clusters arc consistent with the very short timescales 
for cluster dissolution inferred from the previous investi- 
gations of Fall (2004) and Chandar et al. (2005) of local 
SB galaxies 

2. A Model for the Evolution of the gas 

At least two physical parameters are essential to deter- 
mine the evolution of a proto-cluster: the star formation 
efficiency and the timescale of gas dispersal from the clus- 
ter (e.g., Lada & Lada 2003). 

Accurate measurements of the star formation effi- 
ciency, SFE = Mstar/{Mgas + Matar), which require a 
reliable determination of both the gaseous and the stel- 
lar mass content, are in general not available for cluster 
forming regions, nor even in our galaxy. In the latter, the 
inferred SFEs of proto-clusters embedded in giant molec- 
ular clouds range from ^ 10% to 30%. 

The timescale for gas removal, Tg, is even less con- 
strained by empirical data. In SB systems, in particular 
where the rates of SN explosions arc elevated, one should 
expect an efficient heating of the gas followed by its com- 
plete removal from the proto-cluster. The condition for a 
cluster to become gravitationally unbound after a rapid 
gas dispersal (e.g., in a timescale Tg < Tcr) is that the in- 
stantaneous velocity dispersion of the stars at the time of 
gas removal is larger than the escape velocity, or a > v^sc- 
From the equations for cr and Vesc given in Section 1, one 
finds that this condition implies that a system with SFE 
< 50% should become unbound after rapid gas removal, 
(e.g., Lada & Lada 2003). In the opposite case of a slow 
gas removal (r^ < Tcr) even clusters with low SFEs would 
have time enough to adjust adiabatically and expand to 
new states of virial equilibrium and remain bound. 

However, the simple estimates above, do not take into 
account the fact that proto-clusters are actually far more 
complex systems. For instance, the gas is not homoge- 



neously distributed but may be mostly within dense small 
clouds. These clouds will suffer mass loss by ablation, ther- 
mal evaporation and photo-evaporation due to the large 
amount of ionizing photons emitted by the massive stars 
and by the SNe (e.g., Melioli & de Gouveia Dal Pino 2004, 
2005). The diffuse gas will be heated by the SNe, but may 
cool via radiative losses which depend on the square of the 
gas density, so that on one hand the SN remnant (SNR) 
shock fronts may cause the heating and expansion of the 
diffuse gas, but on the other hand, the mass loss from the 
clouds may increase the density of the diffuse gas enough 
to allow an efficient radiative cooling that may delay the 
gas removal. Defining the heating efficiency of the SNe, 
HE, as the fraction of their energy that is effectively stored 
by the gas in the form of enthalpy and is not radiated away, 
then one may expect that an efficient gas expansion and 
removal will occur only when HE becomes close to unity. 

In order to compute the evolution of the gas of a proto- 
cluster embedded in a SB environment considering the dif- 
ferent physical processes mentioned above, we will adopt 
the semi-analytical dynamical model described by Melioli 
& de Gouveia Dal Pino (2004, hereafter MGDP). 

As in MGDP, we assume that an instantaneous burst 
of star formation occurs inside a spherical parent cloud of 
radius R. We adopt a Salpeter initial mass function (IMF) 
with a number of stars with mass greater than 8 M© (i.e., 
the SN number) J\f ~ Q.QlMstari where M^tar is the total 
mass of stars in Mq. 

For simplicity, it is assumed that the expanding SNRs 
and the cold, dense cloud clumps arc uniformly distributed 
in the parent cloud. The clumps have an initial density 
ric ^ 10^, radius Tc ^ 0.02 pc, and mass Mc = 0.015 
M0. Assuming similar characteristics to those observed 
in the Rosette giant cloud in our galaxy, in most of the 
calculations below we will adopt a total mass of clumps 

Mga, ~ 8M,tar- 

The equations that describe the evolution of the gas 
are given in MGDP. In the Section below, we will present 
the evolution of the diffiise gas in a proto-cluster for differ- 
ent values of the SFE until the time when the heating ef- 
ficiency, HE, which is mathematically defined as the ratio 
of the enthalpy flux through the boundary of the system 
to the total SN energy flux: 

HE = 4.i^^ii^;±M^e., (1) 

/■CEsN 

becomes equal to unity. 
3. Results 

Table 1 summarizes the input parameters of the com- 
puted models. The values adopted are compatible with 
those estimated from the observations (e.g., Chandar et 
al. 2005, Fall 2004, Tenorio-Tagle et al. 2005 and refer- 
ences therein). 

Figure 1 shows the time evolution of the ambient tem- 
perature, density, pressure, and the heating efficiency, HE, 
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Fig. 1. Ambient density (n), temperature (T) and pres- 
sure (p) evolution in the left panel for a spherical stellar 
cluster with initial Mstar = 10^ M©, R=10 pc, ambient 
density n=10 cm'^, T = lO"* K, p = 1.38 x 10"" dy cm'^ 
and total gas mass Mgas = 8 x IO^Mq. In the center panel: 
number of clouds {Nc) and clouds radius (rc) evolution; 
and in the right panel: the SN heating efficiency (HE). 



for the proto-cluster of model 1 of Table 1, for which SFE 
= 11%. After 3 Myr, HE rapidly increases to one. Before 
that, the total energy stored in the gas is only a few per- 
cent of the total energy released by the SNe, most of which 
is radiated away. During this period, the mass-loss rate 
from the clouds, mainly due to the process of photoevap- 
oration (see MGDP), keeps the density of the diffuse gas 
high enough allowing an efficient radiative cooling. After 
this time, the clouds evaporated almost completely and 
the supernovae heating efficiency, HE, increased to unity 
in less than 0.3 Myr and this leads to a quick gas expan- 
sion and removal from the system. Note that right after 
3 Myr the gas temperature rises to T ~ 10^ K and thus 
the gas will freely expand through the cluster radius in 
texp = R/vs < 0.1 Myr. These time scales are smaller than 

~ 0.5 Myr so that, once HE reaches the unity value, 
all the gas will be removed from the cluster in less than 
one cluster dynamical crossing time and, since a > v^sc, or 
SFE < 50%, we may conclude that the system will become 
unbound and disperse. 

Figure 2a compares HE for systems with different ini- 
tial radius. As expected, the larger the radius of the proto- 
cluster the longer the system will take to reach a heating 
efficiency equal to unity, but in all cases we observe the 
same trend as in Figure 1, i.e., a rapid change in HE from 
<< 1 to ^ 1 in less than the corresponding Tc,. 0.5, 
1.4, and 2.6 Myr, for the systems with R — 10, 20, and 
30 pc, respectively. Figure 2b shows the evolution of HE 
for different SFEs. All the curves have the same initial 
conditions of Figure 1 except for the initial total gas mass 
(Mgas)- We see that the higher the SFE (i.e., the smaller 
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Fig. 2. (a) Top: HE computed for different initial radius. 
All the curves have the same initial conditions of Model 1 
of Figure 1 (with a SFE = 11%), except for the radius of 
the system. R = 10 pc is represented by the dot-dashed 
lines, R = 20 pc by soHd Hues, and R = 30 pc by dashed 
lines, (b) Bottom: Evolution of HE for different SFEs. AU 
the curves have the same initial conditions of Model 1 
except for the initial gas mass content. SFE = 11% (with 
Mgas = 8 x 10^ Mq) is represented by the dot-dashed 
lines, SFE = 30% (with Mgas = 2.3 x 10^ M©) by sohd 
lines, and SFE = 50% (with Mgas = 10^ Mq) by dashed 
lines. 



the total gas content of the system relative to the total 
amount of mass in stars) more rapidly the system reaches 
a high HE and therefore, earlier is the gas removal. Also, 
since for the investigated range of SFE values, SFE = 10%- 
50%, the timescales for gas removal right after HE reaches 
the unity value are smaller than the clusters dynamical 
crossing times, Tcr ^ 0.5, 0.8, 1 Myr, for the clusters with 
SFE = 11%, 30%, and 50%, respectively, then they all 
should become gravitationally unbound systems. 

Figure 3 adds information to the results of Fig. 2b, 
showing SFE versus the time the that the system takes to 
reach HE = 1 for clusters with different initial conditions 
(Models 1 to 4 of Table 1). the can be also interpreted as 
the time the gas is retained in the proto-cluster possibly 
allowing for new generations of star formation. The curves 
indicate that this timescale in a cluster is very insensitive 
to SFE when the total initial gas mass content is fixed. 
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Fig. 3. SFE versus the time the that the system takes to 
reach HE=1 for clusters with different initial conditions 
(Models 1 to 4 of Table 1). 

The varies less than 30% between SFE = 0.1 and 0.5 in 
each example. Combined with the results of the previous 
figures, Figure 3 also reveals that due to a quick gas re- 
moval right after HE becomes equal to one, typical SB 
proto-clusters are expected to start to disperse after < 5 
Myr. 

4. Conclusions 

We have investigated the evolution of the gas in proto- 
clusters embedded in SB systems and have shown that 
even in the presence of dense, cold clouds that tend to 
inhibit the heating of the interstellar gas by the SNE, the 
timescales during which the gas is retained in these sys- 
tems do not exceed '~ 5 Myr (considering proto-clusters 
with total mass in the range Mstar + Mgas ^ Qy. 10^ — 6 x 
10^ M0 and radii 5 pc < i? < 20 pc). We have also found 
that after this time the gas is removed very rapidly (^ 0.5 
Myr) in timescales smaller than the clusters dynamical 
crossing time and this result is approximately insensitive 
to the initial SFE for a given total gas mass. These very 
short timescales for gas removal combined with values of 
SFE < 50% will make these systems to start to disperse 
after ~5 Myr. 
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